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Abstract 

We discuss the steps to construct Dirac operators which have arbi- 
trary fermion offsets, gauge paths, a general structure in Dirac space 
and satisfy the basic symmetries (gauge symmetry, hermiticity con- 
dition, charge conjugation, hypercubic rotations and reflections) on 
the lattice. We give an extensive set of examples and offer help to 
add further structures. 
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1 Introduction 



Dirac operators which satisfy the Ginsparg- Wilson (GW) relation 

7 5 L>- 1 +£>- 1 7 5 = 7 5 2i?, (1) 

where R is a local operator trivial in Dirac space, define lattice regularized 
theories whose chiral properties are the same as those of the corresponding 
formal continuum theories Jl]^ . The relation in eq. (Q) implies an exact chiral 
symmetry on the lattice [0 . Together with the earlier results on domain wall || 
and overlap |Q fermions, the recent developments provide a rather satisfactory 
understanding of global chiral symmetry on the lattice. There is significant 
progress on chiral gauge theories also Q . 

Two types of solutions to the GW relation are known: the fixed-point (FP) 
Dirac operators of different renormalization group transformations (2) and the 
Dirac operators obtained by using Neuberger's construction ||. In the latter 
case one might even start with the simplest nearest-neighbour Wilson operator. 
It is probable, however that starting with a Dirac operator which satisfies the 
GW relation approximately will reduce the numerical problems with the square 
root in Neuberger's construction and will more than pay back the additional 
overhead^]. In this case and, of course, in any parametrization of the FP operator 
one faces the problem of constructing a Dirac matrix with different fermion 
offsets, gauge paths and a sufficiently rich structure in Dirac space in such a way 
that the basic symmetries (gauge symmetry, 75 hermiticity, charge conjugation, 
hypercubic rotations and reflections) are satisfied. 

This problem was addressed in d = 2 in a recent paper by Gattringer and 
Hip (!(]] . Considering gauge paths with length up to 2 and all 4 matrices of the 
d = 2 Clifford algebra they discussed and solved the symmetry conditions. The 
parameters were fixed then by optimizing for the GW relation. This program 
has been extended to d = 4 in a follow-up paper by Gattringer jflj. The 
first step, i.e. writing down an Ansatz for a generalized Dirac operator which 
satisfies the basic symmetries, is however independent of the method used to fix 
the parameters. 

In this paper we discuss this kinematic problem in d — 4. For any fermion 
offset and for any of the 16 elements of the Clifford algebra we describe the steps 
to find combinations of gauge paths which satisfy all the basic symmetries. We 
construct explicitly paths for all the elements of the Clifford algebra in the 
offsets of the hypercube. We show also how to factorize the sum of paths (i.e. 
writing it as a product of sums) in such a way that the computational problem 
is manageable even if the number of paths is large. 

2 Starting with the Wilson operator the end result after Neuberger's construction is local 
but very broad and has strong cut-off effects jsj]. This gives additional (and related) reasons 
to start with an operator with better properties. 
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We attempted to organize this paper in a way which makes it possible to use 
the results without the need of repeating the straightforward but lengthy algebra 
involved. Care was taken to introduce notations and to fix the conventions. We 
explain with examples how to use the tables and the listed combinations of gauge 
field products in order to build up the matrix elements of the Dirac operator. 

In case the reader wants to consider offsets and paths which do not enter 
our list, we offer an easy-to-use Maple code dirac .maple which is included 
with this paper in the hep-lat archive. The user should specify a Clifford 
algebra element and a gauge path, and the output will help to construct the 
corresponding contribution to the Dirac matrix. 

Let us add a few remarks at this point: 

1) There are many ways to fix the parameters of the Ansatz. As opposed to 
production runs this problem should be treated only once. It is useful to invest 
effort here, since the choice will influence strongly the quality of the results in 
simulations. We plan to optimize the parameters to the FP action. 

2) There are compelling reasons to use the elements of the Clifford algebra be- 
yond 1 and 7 M in the Dirac operator. For an operator D satisfying the GW 
relation in eq. ([!]), Tt(j 5 RD) is the topological charge density |^,[||, i.e the oc 7 s 
part of D is obviously important. Similarly, a cx a^ u term is already required 
by the leading Symanzik condition. 

3) If the parametrization is good, the GW relation is approximately satisfied, 
the operator under the square root in Neuberger's construction is close to 1. An 
expansion is expected to converge very fast in this case. 

4) The basic numerical operation in production runs is Dry, where D is the pa- 
rameterized Dirac matrix and 77 is a vector. The matrix elements of D should be 
precalculated before the iteration starts. Using all the Clifford algebra elements 
and arbitrary gauge paths, the computational cost per offset of the operation 
Dr) is a factor of ~ 4 higher than that of the Wilson action. Keeping the points 
of the hypercube without the outmost corners (65 offsets), say, the cost per 
iteration is increased by a factor of 4 x 65/8 ~ 30 relative to Wilson Dirac 
operator. 

2 Symmetries of the Dirac operator 

We define the basis of the Clifford algebra as Y = 1, 7 M , icr^, 75, 7^,75. We 
use the notation S,V,T,P and A for the scalar, vector, tensor, pseudoscalar 
and axial-vector elements of the Clifford algebra, respectively. Notice that the 
tensor (T) and axial-vector (A) basis elements of the Clifford algebra are anti- 
hermitian. It will later be useful to list the basis elements of the Clifford algebra 
by a single index, as Fa, A = 1, ... ,16. We choose the ordering 

1 , 71 , 72 , 73 , 74 , ^12 , , icri4 , i(7 2 3 , « CT24 , W34 , 75 , 7 l75 , 7275 , 7375 , 7475 . (2) 
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Under charge conjugation, the Dirac matrices transform as 

C7jC- 1 = - 7M . (3) 

Having chosen the T and A basis elements of the Clifford algebra to be 
anti-hermitian, all the basis elements have the property that 

C^T*^- 1 =T . (4) 

We define a sign er by the relation 

75r f 75 = e r r, (5) 

which gives es = ep = £a = 1 and ey = ct = —1. 

Under reflection of the coordinate axis 77, the basis elements of the Clifford 
algebra are transformed as T' = V^"P~ X , where V n = J v j5- (When T is writ- 
ten in terms of Dirac matrices 7^ this amounts simply to reversing the sign 
of 7 I( .) The group of reflections has 2 4 = f6 elements. Under a permutation 
(1,2,3,4) — > (pi,p2,P3,Pi) of the coordinate axes the basis elements of the Clif- 
ford algebra transform by replacing /1 — > p^ and accordingly 75 — > ep 1 p 2 p 3 p 4 75- 

The matrix elements of the Dirac operator D we denote as 

D(n,n';U) a a i, (6) 

where n, a and a refer to the coordinate, colour and Dirac indices, respectively. 
From now on, we suppress the colour and Dirac indices. The gauge configu- 
ration U appearing in D is not necessarily the original one entering the gauge 
action - it could be a smeared configuration. Provided the smearing has ap- 
propriate symmetry properties (as all conventional smearing schemes do) the 
constructions in this paper remain valid, without any extra modifications. 

The Dirac operator should satisfy the following symmetry requirements: 
gauge symmetry 

Under the gauge transformation U^(n) —> U^(n) 9 = g(n) U^in) g(n + fi)^ , where 
g(n) e SU(N) we have 

D(n, n'; U) -» D(n, n'\ U g ) = g(n) D(n, n'\ U) gin') 1 *. (7) 

translation invariance 

Translation symmetry requires that D(n, n + r) depends on n only through the 
n-dependence of the gauge fields. There is no explicit n-dependence beyond 
that. In particular, the coefficients in front of the different paths which enter D 
do not depend on n. 
hermiticity 

D(n,n';U) =j 5 D(n',n,liy j 5 , (8) 
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where f is hermitian conjugation in colour and Dirac space, 
charge conjugation 

D{n,n';U) = C D(n' ,n;U*) T C" 1 , (9) 

where T is the transpose operation in colour and Dirac space. 
It will be useful to combine eqs. (||,|9|) to obtain 

D(n, n'; U) = C7 5 D(n, n'; U*)* 7 5 C~ 1 . (10) 

reflection of the coordinate axis 77 



D{n,n';U) = V^ 1 D{n,n';U v ^>)V v , (11) 

where V v = j^j 5 and n v — n„ if v ^ 77, while n v = —n v . The reflected gauge 
field U Vr > is defined as 

U^(m) = U^m-ftf, (12) 
U^(m) = U v (m), v^r). 

permutation of the coordinate axes 

These are defined in a straightforward way, by permuting the Lorentz indices 
appearing in D. Note that rotations by 90° on a hypercubic lattice can be 
replaced by reflections and permutations of the coordinate axes. 



3 Constructing a matrix D which satisfies the 
symmetries 

To describe a general Dirac operator in compact notations, it is convenient to 
introduce the operator V M of the parallel transport for direction fi 

(%) = UpWSn+w , (13) 
V / nn f 

and analogously for the opposite direction 

(U- M ) = U^(n - A) t <Jn-A,n' ■ (14) 
V / nn' 



Obviously iU^ 




(Note that in terms of these operators the forward 



and backward covariant derivatives are: 9^ = — 1, 9* = 1 — UL and 9*9^ — 
Up + U\ — 2.) The Wilson-Dirac operator at bare mass equal to zero reads: 

£> w = ^E>(^-^) + ^E( 2 -^-^) • (15) 

fj, ij, 
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It is also useful to introduce the operator U(l) of the parallel transport along 
some path I — [li, I2, ■ ■ ■ , h] where U — ±1, . . . , ±4, by 

U(l) = U h U l2 ...U lk . (16) 

In terms of gauge links this is 

(#(0) , = (yh(n)Ui 3 (n + 8 n+run > , (17) 

where n = l\ + . . . + 1^ is the offset corresponding to the path I. Note that 
U(ly = U(l) where I — [—Ik, ■ ■ ■ , —h] is the inverse path. In particular, one has 

(#([2,1,-2])) =U 2 {n)U l {n + 2)U 2 {n+i)U n+l>n , (18) 

for the corresponding staple. 

As another example, the Sheikholeslami-Wohlert (or clover) term jl2j intro- 
duced to cancel the 0(a) artifacts is given (up to a constant factor) by 

«CV UJ([fi, v, -fj,, -v]) + U([v, -fi, -v, fj])+ 

U([-^-v,ii,u]) + U([-v,ii,u,-ii])-h.c.^ . (19) 
We consider a general form of the Dirac operator 

D = ^2r A J2ofu(i). (20) 

A I 

The Dirac indices are carried by = 1,7„, the coordinate and colour 
indices by the operators U(l). The Dirac operator is determined by the set of 
paths I the sum runs over, and the coefficients cf. 

In the case of -Dw, for T = 1 one has I — [1], [— 1], . . . , [—4] and I = [] (the 
empty path corresponding to {/([]) = 1), while for T = 7^: I = \p] and [— /i]. 
In the clover term, eq. (|l9| ) for T = ia^: I — [fx, v, — fi, — . . . (altogether 
6 x 4 = 24 plaquette products). As these well known examples indicate, the 
coefficients cf for related paths differ only in relative signs, which are fixed by 
symmetry requirements. 

Our aim is to give for all T's and offsets r on the hypercube a set of paths 
and to determine the relative sign for paths related to each other by symmetry 
transformations. We give the general rules for arbitrary offsets and paths as 
well. 



Eqs. (ffpC|) imply that the coefficients cf in eq. (20) are real. Further, from 
hcrmiticity in this language it follows that the path I and the opposite path / 
(or equivalently, U(l) and U{ly) should enter in the combination 

r (U(l) + e r U(iy) , (21) 
where the sign er is defined by 7 5 rt 7 5 = e r r, eq. (g). 
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The symmetry transformations formulated in terms of matrix elements in 
the previous section can be translated to the formalism used here. The reflec- 
tions and permutations act on operators U(l) in a straightforward way. Under 
a reflection of the axis r\ one has U(l) — > U(l') where l[ = —li if \k\ = rj and 
unchanged otherwise. Under a permutation (j>i,P2,P3,pa) a component with 
li = ±/i is replaced by li = ±p M , as expected. The number of combined symme- 
try transformations is 16 x 24 = 384. We denote the action of a transformation 
a = 1,... ,384 by T -> r<- a \ U(l) U(l^). Acting on the expression in 
eq. (|2l|) by all 384 elements of the symmetry group and adding the resulting 
operators together, the sum will satisfy the required symmetry conditions for a 
Dirac operator. 

Let us introduce the notation 

d ~( r ' = Jf E r(Q) + erW a) f) ■ (22) 

a 

A general Dirac operator will be a linear combination of such terms. The nor- 
malization factor N will be defined below. The total number of terms in eq. ( p2|) 
is 2 x 384 = 768. Typically, however, the number of different terms which sur- 
vive after the summation is much smaller. It can happen that for a choice of 
starting V and / the sum in eq. ( ^2[ ) is zero. In this case the given path does not 
contribute to the Dirac structure T. 

To fix the convention for the overall sign we single out a definite term in the 
sum of eq. ([22l ) and take its sign to be +1. Denote by Fq, Iq the corresponding 
quantities of this reference term, and by rg = t(Iq) the offset of Iq. This term is 
specified by narrowing down the set {T^ a \ to a single member as follows: 

a) Given an offset r = r(l) = (n, r%, r3, r^) the reflections and permutations 
create all offsets (±r pi , ±r p2 , ±r P3 , ±r p4 ) where (p%,P2 ,P3,P4) is an 
arbitrary permutation. We choose for the reference offset r the one from 
the set {r(/( Q ))} which satisfies the relations roi > ro2 > ^03 > r 04 > 0. 

b) If several T matrices are generated to this offset then choose as To the one 
which comes first in the natural order, eq. (Q). 

c) Consider all the paths {l^ a '} having offset tq and associated with To, i.e. 
r(/( Q )) = ro and = Tq. To single out one path Iq> from this set, we 
associate to a path I = [li, £2, • • • , h] a decimal code d\d,2 ■ ■ - dk with digits 
di = li if li > and di = 9 + k for k < 0. The path with the smallest code 
will be the reference path Iq. (In other words we take the first in lexical 
order defined by the ordering 1,2,3,4, —4, —3, —2, —1.) 

Of course, one can take To, Iq as the starting T and and we shall refer to 
the expression in eq. (|2^ ) as d(To, Iq) to indicate that it is associated to a class 
rather than to a specific (T, I). 
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We turn now to the normalization of d(To,lo). In general, there will be 
K different paths in the set {l^ \ T^ a ' = T ,r(l^) = r }, i.e. corresponding 
to the same offset ro and Dirac structure To. The normalization is fixed by 
requiring that the coefficient of the reference term TqU(Io) is +1/K . 

Consider a simple example explicitly. Let ro = (1,0,0,0), To = i<j\2 and 
l = [2,1,-2]. The starting term, eq. (fl]) is ia 12 (U([2, 1, -2]) - U([2, -1, -2])). 
Applying all the 16 different reflections gives 

8ta 12 (tl([2, 1, -2]) - U([2, -1, -2]) - U([-2, 1, 2]) + U([-2, -1, 2])) . (23) 

Applying all the permutations on this expression results in: 

d(T , l ) = i [ma 12 (U([2, 1, -2]) - U([-2, 1, 2])) + 
16ia 13 (l>([3, 1, -3]) - U([-3, 1, 3])) + 

16ia u (c>([4, 1, -4]) - U{[-4, 1,4]))+...}, (24) 

where only the terms with the offset r — ro are written out explicitly. Their 
total number is 96. The whole generated set has 8 different offsets giving al- 
together 768 terms. Notice the form of the contribution in eq. (|24|). There 
is a common factor (16 in this case) multiplying all the different operators. 
Only the tensor elements of the Clifford algebra enter, since we started with 
a tensor element. Beyond the common factor the path products have a co- 
efficient ±1. These features are general. The number of different paths with 
T = io- 12 and r = (1,0,0,0) is K=2: the paths [2,1,-2] and [-2,1,2]). 
The normalization factor in this case is Af = 16 x 2 = 32, so that one has 
d(T Q ,l Q ) - \ia 12 (U([2, 1, -2]) - L>([-2, 1, 2])) + . . . . 
A general Dirac operator is constructed as 

f(T ,l )d(T ,l ). (25) 

The coefficients /(r , Iq) are real constants or gauge invariant functions of the 
gauge fields, respecting locality, and invariance under the symmetry transfor- 
mations. These are the free, adjustable parameters of the Dirac operator. 

4 Tables for offsets on the hypercube 

Choosing offsets and paths to be included in the Dirac operator is a matter 
of intuition. It is also influenced by considerations on CPU time and memory 
requirements. In Tables |l||| we give examples of reference paths l for offsets 
on the hypercube and general Dirac structure. The first 3 columns give T , lo 
and the number of paths K as defined in the previous section. The 4th column 
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r 


rcf. path Iq 


K 


T's generated 


cx 


1 


D 

[1,2,-1,-2] 




1 

48 


1 


1,0 
1,0 


7i 


[1,2,-1,-2] 


24 


7 1 ,-..,7 4 







[1,2,-1,-2] 


8 


iai2, • • • , ic 3 4 


-1 


75 


[1,2,-1,-2,3,4,-3,- 


-4] 


384 


7 5 


-1/6 




[1,2,3,4, -1, -2, -3, - 


-4] 


384 




-1/6 


7i 75 


[1,2,-1,-2,3,4,-3,- 


-4] 


192 


7i 75, ■ • ■ ,7475 







[2,1,-2,-1,3,4,-3,- 


-4] 


192 









[2,3,4,-3,-2,-4] 




96 








Table 1: Reference paths for different To's for offset (0000). 



r 


ref. path / 


K 


T's generated 


Cx 


l 


[1] 


1 


1 


8, 1 




[2,1,-2] 


6 




8, 1 


7i 


[1] 


1 


7i 


2 




[2,1,-2] 


6 




2 


72 


[1,2,3,-2,-3] 


16 


72,73,74 







[2,1,-2] 


2 


2C712, ZCT13, ZCT14 


4 


2023 


[1,2,3,-2,-3] 


16 


2O23, «cr 2 4, i<T 3i 


-4 


75 


[2,1,-2,3,4,-3,-4] 


96 


75 


4/3 


7i 75 


[2,1,-2,3,4,-3,-4] 


96 


7i 75 





7275 


[1,3,4,-3,-4] 


16 


7275,7375,7475 


-2 



Table 2: Reference paths for different r 's for offset (1000). 



r 


ref. path Iq 


K 


T's generated 


cx 


1 


[1,2] 


2 


1 


24, 6 


71 


[1,2] 


2 


7i,72 


12 


73 


[1,3,2,-3] 


8 


73,74 







[1,2] 


2 


20"12 


—2 


20-13 


[1,3,2,-3] 


4 


2CT13, 2CT14, ZCT23, 2CT 2 4 





20-34 


[1,2,3,4,-3,-4] 


32 


i(T 3 4 


-4 


75 


[1,2,3,4,-3,-4] 


32 


75 


-2 


7175 


[1,2,3,4,-3,-4] 


16 


7175,7275 


4 


7375 


[1,4,2,-4] 


8 


7375,7475 


-4 



Table 3: Reference paths for different r 's for offset (1100). 
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1 


VAT nof h / „ 

IcI. pcLlll to 




1 fa gCIlCIcLlCU 


cx 


1 


[1,2,3] 


6 


1 


32, 12 


7i 


[1,2,3] 


4 


71,72,73 


24 


74 


[1,2,4,3,-4] 


24 


74 







[1,2,3] 


4 




-8 


i(Ti4 


[1,4,2,-4,3] 


8 


i(Ti4, ZCT24, «CT34 





75 


[1,4,2,-4,3] 


12 


75 


-8/3 


7i75 


[1,4,2,-4,3] 


8 


7175,7275,7375 


8 


7475 


[1,2,3] 


6 


7475 


-4/3 



Table 4: Reference paths for different Tq's for offset (1110). 



r 


ref. path lo 


K 


T's generated 


cx 


l 


[1,2,3,4] 


24 


1 


16, 8 


7i 


[1,2,3,4] 


12 


7 1 ,-..,7 4 


16 


20-12 


[1,2,3,4] 


8 


ia 12 , . . . , icr 3 4 


-8 


75 


[1,2,3,4] 


24 


7 5 


-2/3 


7i75 


[1,2,3,4] 


12 


7175, . . . ,7475 


8/3 



Table 5: Reference paths for different Tq's for offset (1111). 



gives those Clifford basis elements which are generated in eq. ( p2| ) to the offset 
ro- 

It is also of interest how the given Dirac operator behaves for smooth gauge 
fields, i.e. to obtain the leading terms in the naive continuum limit. The 
expression in cq. ( p2|) generated by given Tq and Iq contributes in this limit to 
one of the expressions (of type S,V,T,P,A) below 

(c s + c s d 2 ) ; cyjf.df,; or-a^F^ ; 

cv^ihe^ p <jF^F pa ; CAH^5^^ P ad u F pcr . (26) 

Here is the covariant derivative in the continuum. The coefficients cx(ro, lo), 
which determine the continuum behaviour of the Dirac operator, are presented 
in the last column of Tables For To = 1 the first and second entries 

correspond to cs and cs, respectively. We list their meaning below. 
Introduce the notation 

C X = X)/(r 0l Jo)e*(To,lo). (27) 
The bare mass mo is given by 

C s = m . (28) 
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The normalization condition on the D ~ 7^9^ term in D gives 



CV 



1 . 



(29) 



The 0(a) tree level Symanzik condition reads 



C s + C T = . 



(30) 



The coefficient cp is interesting if the parametrization attempts to describe 
(approximately) a GW fermion. In this case it is related to the topological 
charge density q(n) = Tr(7 5 Di?)„ n , 



(If /(ro,^o) is not a constant, then its continuum limit is understood here.) Of 
course, in the expressions above Tq should be of the corresponding type (S,T, 
. . . ). The quantity R e (ro) is given by R(n, n + tq; U — 1) of the GW relation 



5 The list of factorized contributions 

As the tables |l|-^| show the number of paths, in particularly for the offset tq = 
(0, 0, 0, 0), is large. It is important to calculate the path products efficiently. An 
obvious method is to factorize the paths, i.e. to write the sum of a large number 
of paths as a product of sums over shorter paths. We shall try to factorize in 
such a way that these shorter paths are mainly plaquette, or staple products. 

Define the plaquette products (i.e. the operators of parallel transport along 
a plaquette) in the following way: 



where = ±1, . . . , ±4. Their hermitian conjugate is given by (Pi lt i 2 ) — Pl 2 ,h- 
Reflections and permutations act on them in an obvious way. 
We also define the staple products as 




(31) 



eq-|). 



Pi u l 2 =U([h,h,-h,-h]), 



(32) 



U([k,l 2 ,-h]). 



(33) 



We have (Si u i 2 y = Si x -i 2 . 



To describe the shortest path to an offset we introduce the notation 



V h = U([h]), 

Vl u l 2 = U([h,l2]), 
Vl u l 2< l 3 =U([h,l2,k]), 

Vh,l2,l3,h — ^([^1^2, h]) ■ 



(34) 
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Introduce the following linear combinations transforming in a simple way 
under reflections: 



p(++) — p ip ip ip 

Paw ^ P^-. — v ~f Mi w — M) — w ' (35) 

p( — H) — P ip _ P _ P 

/^f * fJ,,V ^ - 1 fJ,, — U x — x — /X,— f ) 

p( ) — P _ P _ P IP 

The signs in the superscript denote the parity for reflections of the /x, i/ axes, 
respectively. Hermitian conjugation acts as interchanging both upper and lower 
indices, permutations simply by /i — * p M , ^ — > p„. 

It is also useful to denote combinations which are symmetric/ antisymmetric 
with respect to interchanging the axes: 

P { »7 m) = P^ +) + P& +) = Pfr +) + h - c - . (36) 

p(as) _ p(--) _ p( ) _ p( ) _ Vl c 

jLiV " 



For the staples we write 



Cf(w+) _ o , o 

^ — f_L k', - M * u — ^ 

o(w-) _ r, Q 
^ — ^ IS] — /J' ^ — 

The subscript ±/k denotes the direction of the staple, the superscript specifies 
the plane /xj/ and the parity in v. 
For a line of length 2 we have 

^ ( + +) = + V»,-» + ^ + V-^-v , 

= V^ v - V^-u + V-w - V-^-u , (38) 
vfc +) = v MtV + V M ,- V ~ V-^. v - V-^-u , 



and analogously for longer lines. In analogy to eq. ( |3q ) we introduce the com- 
pletely (anti)symmetric combinations V^T.. , V^}. for the line products. 
Define also the bent rectangle, 

BhMa — U([h,l2,h, —h, —h, ~h]) , (39) 

and Bvp a "* which is odd in all its arguments (as a sum of 8 terms, defined 
analogously to Pji v ' ) . Similarly, we introduce the "4d plaquette" 

Qh,h,lsM = U([h,l2,h,U, —h, —fo, —h, —U}) , (40) 
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and the odd combination q\ , , ) . 

To simplify the expressions below, we require that all the directions entering 
Pt lt l 2 > Sh,hi Vj 1; i 2j ..., etc. are different, i.e. they are taken to be zero if e.g. 
\h\ = U ' 

Below we list d(T ,l ) for different choices of T and l . The notation ^ 
indicates that all indices in the corresponding sum are taken to be different. 

5.1 The offset r = (0, 0, 0, 0) 

r = Mo = [] 

1. (41) 

r = Mo = [i,2, -i,-2] 

^ E P £r ] • (42) 

r = 7i, lo = [1,2,-1,-2] 

^E'^(^ +) -h.c) . (43) 

r = i<7i2, = [1,2,-1,-2] 

^E^ P ^ S) - ( 44 ) 
r = 75, lo = [1,2,-1,-2,3,4,-3,-4] 

Jjg^ 5 ^ -^nvpaP^v' ] Ppa^ ■ (45) 

r = 75, Zo = [1,2,3,4,-1,-2,-3,-4] 

2g4^ 5 e V>vp<7Q\ lvp(y ^ • (46) 

\xvpu 

r = 7175, - [1, 2, -1, -2, 3, 4, -3, -4] 
1 ^ ' 1 



192 E ' VT5 2 W (^"^ + h.c.) . (47) 
r - 7175, = [2, 1, -2, -1, 3, 4, -3, -4] 

E 'v^W ( P i; +)P ^ S) +h-c) . (48) 

[iv pa 

r = 7i75, = [2,3,4,-3,-2,-4] 

| E ' VT5 W (s<~ "> + h.c.) . (49) 



96 

pvpa 



12 



5.2 The offset r = (1,0,0,0) 
r = i, lo = [i] 



E + v -^- ( 5 °) 



r = 1, lo = [2, 1, -2] 



6 ^ 

fJ,U 



+ S^) . (51) 



r = 7i, io = [i] 



E^(^-^-f)- (52) 



r = 7 i, Jo = [2, 1,-2] 



6 



1 v^' 
16 



r = 72, lo = [1,2,3,-2,-3] 

E' > (Vi; +) - h.c.) . (54) 
r = i<7i2, «o - [2,1,-2] 

^EV^^-^-H • (55) 

r = *CT23, = [1,2,3,-2,-3] 

^ E'^4 (^ s) + p ^ v p - h - c -) • ( 56 ) 



E' ' ' 

r =75, /o = [2,1,-2,3,4,-3,-4] 



1 ' 1 

96 75 2 



w (^' _) ^ S) + + h - c -) • ( 5? ) 
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r = 7175, = [2,1,-2,3,4,-3,-4] 

^ E ' 7m75 ^ W (sjr^pM P^S^ + hx.) . (58) 

\xvpo 

r = 7275, Z = [1,3,4,-3,-4] 

^ E ' 7^75 ^ W + V, + h.c.) . (59) 

5.3 The offset r = (1,1,0,0) 

r = i,/ = [i,2] 

^E^ ym) - (6°) 

r =7i, «o = [1,2] 

^E'7.(^ +) -h.c) • (61) 

r =73, i = [1,3,2,-3] 

^ E'7p + V-pSl"-) S^V V S^V- V h.c.) . (62) 

fu/p 

T = ia 12 , l = [1,2] 

\ E • (63) 

r = i(7i 3 , 2 = [1,3,2,-3] 

J E ' ia w tys^ - ^-^^ _) + - l-^^ _) - h.c.) . (64) 

r = i<734, = [1,2,3,4,-3,-4] 

^ E'-p4(^ ++)p ^ s) - h - c -) ■ (65) 
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r = 75, lo = [1,2,3,4,-3,-4] 



^75 E'^W (^--)p^)+h.c.) . 
r =7175, Jo - [1,2,3,4,-3,-4] 

^ E (Vfi-)P<£> +h.c.) . 

r - 7375, Zo = [1,4,2,-4] 



5.4 The offset r = (1,1,1,0) 
r = l,/ = [1,2,3] 



r = 7i, Jo = [1,2,3] 



I \ " T/(sym) 



6 

H<v<p 



iE^(^ +) + ^- } ) 

P^P 

r = 74, Jo = [1,2,4,3,-4] 



r = i<7i2, Jo = [1,2,3] 

T = i<Tu, Jo - [1,4,2,-4,3] 
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To = 75, lo = [1,4,2,-4,3] 



(74) 

To = 7i75, = [1,4,2,-4,3] 

-VnSf^V-a - V- a S^-% - V-nSf'-lV-r - V-rS^V-n + h.C.) . 

(75) 

To = 7475, h = [1,2,3] 

\ E ' 7*75^ W (y^ P ~~ ] + h.c) . (76) 

5.5 The offset r = (1,1,1,1) 

r = l,/ = [1,2,3,4] 



^ ( 23 y 4 m) ■ (77) 



To = 7i, Zo = [1,2,3,4] 

^eV(^u ++) -H- (78) 



\xvpo 

r = i<ri2, Jo = [1,2,3,4] 

1 v — > ' 



£ ^(^ ++) -h-c.) . (79) 
To = 75, = [1,2,3,4] 

^75 E ' ^ W ( V ^p°~ ] + h - c -) • ( 8 °) 



\xvpu 



To = 7i75, i = [1,2,3,4] 



12 ^ 

livpa 



7^75 Wa + h.c.) . (81) 
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6 Concluding remarks 



We have given a construction of a generalized Dirac operator with the required 
symmetry properties. For every offset on the hypercube and every element of 
the Clifford algebra we have provided at least one example of a contribution to 
the Dirac operator. 

Theoretical or computational considerations should determine the impor- 
tance of a particular term, which could require the reduction or extension of 
the set contained in Tables 1-5. We plan to optimize the parameterization to 
approximate the fixed-point Dirac operator, in order to have good scaling and 
chiral properties. 

Acknowledgment: We thank Philipp Rufenacht for useful discussions. 
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